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Abstract 

Rankine-Hugoniot shock adiabats are calculated in the pressure 
range 1 Mbar-10 Gbar with two atomic-structure models: the atom in 
a spherical cell and the atom in a jellium of charges. These quantum 
self-consistent-field models include shell effects, which have a strong 
impact on pressure and shock velocity along the shock adiabat. Com- 
parisons with experimental data are presented and quantum effects are 
interpreted in terms of electronic specific heat. A simple analytical es- 
timate for the maximum compression is proposed, depending on initial 
density, atomic weight and atomic number. 

1 Introduction 

The study of radiative properties of high-energy-density plasmas has many 
applications in inertial confinement fusion and astrophysics. It requires 
knowledge of the plasma properties in extreme conditions of density and 
temperature. For example, the pressure in the center of Jupiter is believed 
to reach 100 Mbar and temperature 25 kK, while in the center of white 
dwarfs pressure can exceed Gbars at temperatures of 10^ K. In the labo- 
ratory, laser-based shock-wave experiments attain pressures of hundreds of 
megabars, obtained via the absorption of an intense laser pulse. The thermo- 
dynamics and the hydrodynamics of hot dense plasmas can not be predicted 
without a knowledge of the equation of state (EOS) that describes how a 
material reacts to pressure. For instance, the theory of stellar evolution is 
affected by uncertainties in EOS. After being predicted, brown dwarfs have 
recently been observed. Their internal structure and cooling time depend on 
the details of the EOS at densities approaching solid density and at temper- 
atures of a few eV, conditions obtainable in laser experiments pQ. Therefore, 
the need for suitable EOS of high-energy-density matter becomes crucial. 
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A material subjected to a strong shock wave is compressed, heated and 
ionized. As the strength of the shock is varied for a fixed initial state, 
the ensemble of the pressure-density final states of the material behind the 
shock, named shock adiabat or Rankine-Hugoniot (RH) curve, depends on 
the EOS of the matter, which must be determined from theory. At inter- 
mediate shock pressures, when the material becomes partially ionized, the 
EOS depends on the precise quantum-mechanical state of the matter, i.e. 
on the electronic shell structure. Therefore, for the past twenty years there 
is a great interest in the physics of bound levels in high-energy-density plas- 
mas and quantum self-consistent-field (QSCF) models [2] are replacing the 
older Thomas-Fermi (TF) approach [3]. We present RH curves calculated 
using our QSCF-based EOS code, ESODE (Equation of State with Orbital 
Description of Electrons). The ionic contribution to the EOS is described 
by an ideal gas with OOP (One Component Plasma) corrections, and the 
cold curve (T=0 K isotherm) has been obtained in most of the cases from 
Augmented Plane Wave (APW) calculations [3]. 

ESODE has two treatments for the thermal electronic contribution to 
the EOS: Average atom in a Spherical Cell (ASC) or Average atom in 
a Jellium of Charges (AJC). In both cases, bound electrons are treated 
quantum-mechanically. In the ASC model, all electrons are confined within 
a Wigner-Seitz sphere and free electrons are described in the TF approxima- 
tion. In the AJC model, bound and free-electron wavefunctions can extend 
outside the sphere, where the plasma is represented by a uniform electronic 
density (jellium) neutralized by a continuous background of positive charges, 
representing ions. The QSCF models (ASC and AJC) are described in Sec. 
[2j Shock adiabats calculated from those models are presented, analyzed and 
compared to traditional TF model and to published experimental data in 
Sec. [3l General features of RH curves, like the interpretation of shell effects 
in terms of specific heat, are discussed in Sec. E] as well. An analytical 
estimate for the maximum compression rate is proposed in Sec. HI Finally, 
the dependence of shock velocity on particle velocity is analyzed in Sec. O 

2 The quantum self-consistent-field models 

In the present work, we consider the regime identified as strongly coupled 
(non-ideal) plasmas, characterized by a high density and/or a low temper- 
ature. In such plasmas, ions are strongly correlated, electrons are partially 
degenerate, and the coupling parameter F, ratio of Coulomb potential energy 
and thermal energy, is greater than unity. The main task is to evaluate the 
thermal and compressional excitations of the electrons. The determination 
of the average electronic charge density relies usually on the Density Func- 
tional Theory (DFT), of which a well-known example is the Thomas- Fermi 
(TF) model It contains essential features to characterize the material 
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properties in extreme conditions and is expected to be most reliable at ex- 
treme conditions when the detailed influence of the electronic structure does 
not play any role. Despite the omission of quantum structure, TF model 
typically predicts electronic densities in broad agreement with detailed ap- 
proaches. However, as the pressure or temperature is raised, successive shells 
of electrons are delocalized and the effect of such phenomena is represented 
only in an averaged manner by the TF model. 

2.1 Atom in a Spherical Cell (ASC) 

Atoms in a plasma can be represented by an average atom confined in a 
Wigner-Seitz (WS) sphere, whose radius r^js is related to matter density. 
Inside the sphere, the electron density has the following form: 

n{r) = V',(f)' + ^^5^^^[Fi/2(-F(r),/2-F(r),x) 

b 

+xF^,2{-V{r),li-V{T),x)l (1) 

where 



Ji{x,y) = — FT w, ^1 and F„/2 a,x,cj = / — — —dy 

l + exp[{x-y)/kBT\ ' 1 + exp[y - x] 

(2) 

are, respectively, the Fermi-Dirac population and the modified Fermi 
function of order n/2. The first term in Eq. ([T]) corresponds to the contri- 
bution of bound electrons to the charge density, while the second term is the 
free electron contribution, written in its semi-classical TF form. The quan- 
tity eh is the energy of a bound orbital and ^pb the associated wavefunction 
calculated in the Pauli approximation [5] , in which only first-order relativis- 
tic corrections to the Schrodinger equation are retained. We introduce the 
notations x = ^bT/Eq, Eq being the rest-mass energy of the electron and 
V{r) = V{r)/{kBT), where V{r) is the self-consistent potential: 

V{r) = --+ -^d\' + y.,[n], (3) 
Jo \r — r'\ 

Vxc being the exchange-correlation contribution, evaluated in the local 
density approximation [6j. Last, the chemical potential /i is obtained from 
the neutrality of the ion sphere: 

IS 

n{r)AiTr'^dr = Z, (4) 







and /i = fi/{kBT). Equations ([I]), ([3]) and @ must be solved self- 
consistently provided that bound orbitals are obtained from the Schrodinger 
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equation at each step. The electronic pressure Pe [7] consists of three con- 
tributions, Pe = Pb + Pf + Pxc, where the bound-electron pressure Pf, is 
evaluated using the stress-tensor formula 



(5) 

Uh representing the radial part of the wavefunction ipi, multiplied by r. 
The free-electron pressure Pf reads 

Pf = ^{kBTf/^[F3/2{-V{r^s),fl-V{r^s),x) 

+^F,/^{-V{r^s),fl-V{r^s),x)] (6) 

and Pxc is the exchange-correlation pressure evaluated in the local den- 
sity approximation |6j. The choice of the boundary conditions plays a major 
role in the expression of pressure due to the fact that the energy of an orbital 
depends on the value of the corresponding wavefunction at the boundary of 
the WS sphere. In our model, the wavefunction behaves like a decreasing 
exponential at the boundary. The internal energy in the AJC model is 

Ee = y^m^-- n{r) —^d^rd\'+E^c- n{r)Vxc{n{r))d^r, 

■ ^ JO Jo \r — r'\ Jo 

(7) 

where E^c is the exchange-correlation internal energy and qi the popu- 
lation of state i (either bound or free). The first term in Eq. ([7]) can be 
expressed by 

i b 

Ef^p is the potential energy 

Ef,p = ^(^^^)'^' r\F,/,{-V{r),fl- V{r),x) 
^ Jo 

+xF3/2{-Vir),fl - V{r), x)]V{r)d'r (9) 
and Ef^if. the kinetic energy 

Ef,k = -^^^^p'" r\Fy,{-V{r),fi- V{r),x) 
^ Jo 



efe , 1 + 1 + f 



+xF5/2{-V{r),fL-V{r),xWr. (10) 
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Here only non-relativistic calculations are performed, which correspond 
to E'o ^ in Eqs. ([1]), ([5]) and ([6]), to make comparisons with the non- 
relativistic model described in Sec. 12.21 Equations ([T]) to (llOh enable one to 
include relativistic effects without solving Dirac equation. 

2.2 Atom in a Jellium of Charges (AJC) 

In order to go beyond the TF approximate treatment of continuum electron 
charge density, it is necessary to use a full quantum-mechanical description 
of the continuum states and to consider that both bound and free orbitals 
can extend outside the WS sphere, which requires the definition of the en- 
vironment beyond the central average ion. One way to address this is the 
jellium model (or electron-gas model), i.e. a uniform electron density (— n_(_) 
neutralized by a positive background n+ simulating the ionic charges. The 
model relies on a method proposed by Friedel [U [9] to treat the electronic 
structure of an impurity, represented by a spherical potential of finite range, 
in an electron gas. It has been further developed by Dagens [10^ [TT] and 
Perrot [121 [13l . The problem is reduced to the response of the electronic 
density to the immersion of a positive charge Z into the jellium. These 
considerations lead to the following form of the electron density: 

n{r) = y2Meb,fJ,)M^ +E/ M^, i^eA^) de (11) 
b I 

with 

V{r) = + / -^dV + V,,{r) - V^,{-n+), (12) 

^ Jo \r — r'\ 

where Roo » fa, Va being the radius of the cavity. The potential V{r) 
is determined, as in the ASC model, in a self-consistent way. The ionic 
density is modeled by 

n+(r) = |° ^^^^"^ (13) 
[ Uj^ for r > Ta- 

The expression of internal energy is [14j 

E[n,n+] = Z*[ek + e^c]{-n+) + ^E[n,n+\- KrPiVa (14) 

with 



4 o . d\nZ* 



T 



3 51nT 



Va,T 



^ r Roc 

and Pi = -t / 47rrV(r)dr. (15) 

Jr„ 
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Quantities and e^c are the kinetic and exchange-correlation energies 
per free electron and A£^[n,n_|_] is the immersion energy, i.e. the energy 
change resulting from the immersion of an ion in the jellium: 



Roc fRoo 7 



A£'[n,n+] = / {ek[n{r)] - ek[-n+\)d r - I —{n{r)+n+{r))dr 



JO 



r 
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/o JO 

Raa 

'3, 



+ / (exc(n(r)) - e^c(-n+))(i'^r. (16) 

JO 

The pressure is obtained by 



Pe = [^fc + -Pxc](-n+)+n+[/(ra) + (l-A„jPi, where A„„ 



(91nw 



(17) 



and 



Z 1 



U{r) = h-/ 47rr'^(n(r') +n+(r'))dr' + / 47rr'(n(r') + n+(r'))dr'. 

r r Jq Jj. 

(18) 

The term P^, is the pressure of a free-electron gas. The pressure in Eq. 
()17p is rigorously the derivative of energy with respect to volume. This is 
the main difference with Liberman's model [151 [2], The major difficulty with 
these models is that the average ionization is not well defined when the outer 
electrons are delocalized. The only way to make the formalism variational 
is to specify the ionization in the AJC model. In other words, the ques- 
tion is how to define the residual electron density (— ri+) far away from the 
point where the positive charge is introduced into the jellium. A convenient 
choice is |14j Z*{v) = Z^p{v), Z^p being the TF ionization. Then deriva- 
tives of ionization with respect to volume and temperature can be obtained 
analytically, using the numerical fit proposed by More jT6|. At first sight, 
ASC and AJC models seem to be very different concerning the modeling of 
the environment of the atom, isolated and confined in the ASC model, and 
immersed in an infinite effective medium in the AJC model. However, when 
the electronic structure is calculated in the TF approximation, these models 
are equivalent if the jellium density (— n+) is equal to TF density evalu- 
ated at the WS radius. Such a property comes from a variational principle 
(minimization of the free energy). 

2.3 Ionic contribution and cold curve 

The adiabatic approximation is used to separate the thermodynamic func- 
tions into electronic and ionic components. The total pressure can be writ- 



6 



ten: 



P{p,T) = P,{p) + Pi{p,T) + Pt{p,T), (19) 

where Pt{p,T) = Pe{p,T) — Pe{p,0), is the thermal electronic contribu- 
tion to the EOS. The quantity Pi is the ionic pressure and subscript "c" 
characterizes the cold curve obtained from APW [4J simulations or using 
the Vinet [T7] universal EOS. Equation p9|) also holds for internal energy. 
In order to take into account non-ideality corrections to the thermal motion 
of ions we use an approximation [18] based on the calculation of the EOS of 
a One-Component Plasma (OOP) by the Monte Carlo method [19]. The ion 
contribution can be obtained using a formula based on the Virial theorem: 

Piip, T) = pkeT + |AEi(p, T), (20) 

where 

^E,{p, T) = ksT [r3/2 ^ j^-f^, - air] (21) 

sindEi{p,T) = 3kBT/2+AEi{p,T) with ai = -0.895929, aa = 0.11340656, 
ag = -0.90872827, 04 = 0.11614773, h = 4.666486, 62 = 13.675411, 
63 = 1.8905603 and 64 = 1.0277554. Such corrections lead to small dif- 
ferences in the RH curves as shown in Sec. \E\ 

3 Rankine-Hugoniot shock adiabats 

3.1 Definitions and numerical investigations 

The initial state of the plasma is characterized by a density po, a temperature 
To, a pressure Pq, and an internal energy Eq. D is the shock velocity, u 
the matter velocity, P and E are respectively the pressure and internal 
energy behind the shock front. These variables obey the following Rankine- 
Hugoniot relation [20] : 

hp + P,)(^-^)=E-Eo, (22) 
2 \po P^ 

The RH curve can be obtained by solving Eq. ()22p at each temperature 
step. Calculations have been performed for T < 6.4 keV for beryllium 
(Be, Z=4), aluminum (Al, Z=13), iron (Fe, Z=26) and copper (Cu, Z=29) 
calculated from pure Thomas-Fermi EOS, ASC model, and AJC model. The 
cold curve has been calculated using Vinet universal EOS [TT] for Be and 
APW [U method for Al, Fe and Cu. Results from the QSCF models differ 
strongly from the TF approximation and slightly from eachother. In the 
case of Al, which we use as an example (see Fig. [1]), the difference between 
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the theories appears for P > 3 Mbar and AJC gives higher pressures than 
the other models for 2< rj < 3.5, where r] = p/po is the compression rate. 
For 7] > 3.5, all QSCF models give lower pressures than TF model. 

3.2 Quantum orbital effects 

Our models emphasize the thermodynamic domain where RH curves strongly 
depend on the electronic structure, i.e. beyond four times solid density 
where the shoulders (double in the case of Al, Fig. [1]) correspond to ion- 
ization of successive shells. These shoulders can be explained by the com- 
petition between the release of internal energy stored in the shells and the 
free-electron pressure. When ionization begins, the energy of the shock de- 
populates the relevant shells and the material is very compressive. However, 
the pressure due to the free electrons eventually dominates and the mate- 
rial becomes less compressive. Both models show compression maxima in 
the range 5po — 6po. In this region, the electrons from the ionic cores are 
being ionized. The shock density increases beyond the infinite pressure of 
4po in the electron ionization region. As ionization is completed, the plasma 
approaches an ideal gas of nuclei and electrons and the density approaches 
the fourfold density ApQ. For Be, all the models yield a single density max- 
imum, corresponding to the ionization of the K electron shell. For Al (Fig. 
[1]), there are two density maxima corresponding to the K and L electron 
shells, each followed by density decreases. For Fe and Cu there are den- 
sity maxima or inflexions corresponding to the K, L and M electron shells. 
The L-shell ionization feature gives the largest density increase. In the case 
of Cu, the ASC model exhibits a kind of discontinuity around 40 g/cm^, 
due to pressure ionization of 2p orbital. Such a sharp increase of pressure 
does not exist in the AJC model, because this model is consistent from the 
thermodynamic point of view. This is due to the fact that the pressure is 
rigorously obtained as a derivative of the free energy, and shape resonances 
are carefully taken into account in the quantum treatment of free electrons. 
Such features lead to a continuous disappearing of a bound state into the 
continuum, which is a suitable treatment of pressure ionization. In fact, the 
non-monotonic character of thermodynamic variables stems from the eigen- 
energies of the orbitals themselves, which exhibit the oscillations as well. 
The first density for which decompression occurs is named "turnaround" 
point. The pressure differences from the quantum mechanical theory can be 
explained by examining the electronic heat capacity per particle predicted 
by the two theories along the RH path. At low temperature, the electronic 
heat capacity depends on the number of electrons that can be excited around 
the Fermi energy and the TF theory predicts a smooth increase since the 
density of states in this model is a monotonic function of energy. Therefore, 
the electronic specific heat per particle 
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CI' 



where C' 



lel 



d[E{p,T) - E,{p,T)] 
dT 



(23) 



3fcij/2 



V 



displays the signature of the ionization of successive shells. Figure [2] 
represents C*^' along the RH curve in the TF and in the ASC models for 
Al. Both theories show the effect of the coulomb attractive potential of the 
nucleus binding the electrons, represented by the peak around 300 eV for TF 
theory and 100 eV for ASC model. After the first reduction in density on 
the RH curve, there are 11 free electrons and 2 bound electrons remaining 
in the Is orbital (K shell), which is far away from the energy zero (1.5 keV 
at 100 eV). As long as temperature is not sufficient to ionize those two 
electrons, the specific heat tends to an asymptote corresponding to an ideal 
gas of 11 independent particles. When both Is electrons are ionized (after 
a "threshold" temperature), there is a sudden break in the specific heat, 
which tends to an ideal gas of 13 independent particles. This phenomenon 
is a kind of "Schottky anomaly", i.e. an enhancement in the specific heat 
(see Fig. [5]). This effect is not as important for the 2s and 2p bound states 
(L shell), since their energy levels are not as far from the continuum (a few 
tenth of eV). We note that the same phenomenon occurs for Fe; in that 
case, after two "Schottky anomalies" , the electronic part of the specific heat 
tends to an ideal gas of 26 electrons. 

3.3 Comparisons with experimental data 

Experimental data concerning Be, Al, Fe and Cu [21] have been collected for 
comparison with the EOS models presented. Several experimental methods 
have been used to generate well-defined shock states: gas guns for pressures 
up to 5 Mbar, explosive-driven spherical implosions, laser-driven plane waves 
for generating shocks up to 10 Mbar, and through underground nuclear ex- 
plosions where pressures of a few Gbar can be attained. Superimposing 
these data should yield a single smooth RH curve. The maximum pressures 
reached in the experiments are: 18 Mbar for Be, 4000 Mbar for Al, 191 Mbar 
for Fe and 204 Mbar for Cu. It is clear that there is difficulty discriminating 
between models, since there are very few available data for the region of 
interest, i.e., above 100 Mbar). The error bars associated with the highest 
pressures for Al are too large to provide insight into the existence of shell 
effects. Analysis of the computational results shows that the deviation from 
the experimental points of Al can not be explained only by shell effects. In 
the region where most of the experimental points are available, the role of 
the cold component is important, while shell effects begin to play a signif- 
icant role after the matter is compressed, which occurs near the limiting 
compression r/ ~ 4, and heating begins. For Fe and Cu, results from AJC 
model are found to be in better agreement with experimental points than 
the hybrid model ASC. Finally, we note that it is very difficult to relate gas 



9 



gun measurements with a discussion of shell effects as the pressures are too 
low to exhibit these effect. 



4 Evaluation of maximum compression 

The maximum compression rate depends on the choice of the boundary con- 
ditions of the wavefunctions described in Sec. 12. 1[ Furthermore, it depends 
on the calculation of the ionic part: ideal gas (IG) with or without OOP 
corrections (Eqs. (|20ti2ip ). For instance, the maximum rate for Al in the 
ASC model with IG is 4.901 and with IG and OCP corrections 4.931. OCP 
corrections systematically increase the maximum compression rate. As can 
be checked in Fig. [T]the maximum compression attainable by a single shock 
is larger than 4 and occurs at finite pressure. This phenomenon is due to 
the draining of internal energy in internal degrees of freedom, via ionization. 
An analytical expression for the maximum compression attainable by a sin- 
gle shock in any material from any initial state, except those with gaseous 
densities, can be formulated. The total energy can be written as the sum of 
kinetic and potential energies and Ep. Neglecting exchange-correlation 
contribution, the virial theorem enables one to relate pressure, kinetic en- 
ergy and potential energy 3P/p = 2Ek + Ep. The compression rate rj = p/ pQ 
for the standard RH curve, Pq = 0, po solid density and To=300 K, can be 
written 



r/ = 4 + ^ with 5Ek = Ek - E^^ and 6Ep = Ep - Ep^ . (24) 

^ + ^m; 

At high compression, assuming that Ek » E^^ and that all the electrons 
have been ionized and have a kinetic energy equal to the Fermi energy, one 
can write: 

6Ek = z\i3n'Zp^f/'al (25) 

where ao=52. 9177208319 10^^'^ cm is the Bohr radius, p is in g/cm^, 
and A in g. Next, at high compression the excess potential energy can be 
estimated as the Coulomb interaction energy of two ionic spheres 



1 Z' 



2 



3A 



A-KNAaf 



1/3 



(26) 



SEn = and r,, 

where r^g is the WS radius equal to the cavity radius and we set 
47reo = 1. Therefore, using Eqs. (|25|) and (|26|) . which are relevant for 
a strongly coupled gas of degenerate electrons, the maximum compression 
rate rjm obeys 
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with 



l + C(/Oo,^,^)r/m^ 



C(/5o,^,^) = 37r(2iV^)i/3ao(|L)'/'. (28) 



The solution is: 



1 - (-i)^(c-c^)2cVMC)^i _ _ 32C3 - 1 .3 



(29) 

where Cc = 0.314980262473, d is Heaviside function, 

1 2i°/3 AV3f/-) 



and A(C) = -7 + IGC^ + \/49 + 1824(3 + 256(6. Figure El confirms 
the fact that the maximum compression is always smaller than 7 [52], and 
strongly dependent on the density po- Note that Eq. (f29]) can not be applied 
for gaseous ambient elements. Neglecting cohesive and dissociation energies, 
and using fits for the total ionization energies, Johnson ^22j has proposed 
an analytical formula for the maximum compression rate. It seems difficult, 
however, to discriminate between his approach and ours with the values dis- 
played in Fig. [3] and calculated with our QSCF models, since the disparity 
of values is as large as the difference between the two models. However, 
the present calculation does not rely on a particular form of the EOS, as 
in Ref. [22j. On the contrary, here the maximum compression predicted by 
our model is higher for Fe than for Al, which is consistent with the results 
presented in Ref. However, in Ref. |2j the maximum compression seems 
to increase with Z; according to Johnson's model and ours, this is not true 
stricto sensu but is only a global trend. 

5 Shock and particle velocities 

The particle velocity and shock velocity are: 



V pop \IPo{p-po) 

Expressions in Eq. (j3ip are generic as they do not involve, a priori, any 
explicit relation D{u). For metals, it is often assumed that the relation is 
quasi- linear. However, we find that the (D,u) relationship is almost linear 
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over a wide range of densities, except near n = where we could not perform 
the calculation. It is easy to check from expressions in Eq. (j3ip that the slope 
should be ~ 4/3, corresponding to the ideal gas. However, when looking 
carefully at the first and second derivatives of shock velocity versus particle 
velocity, one finds that the behaviour of shock velocity is more complicated, 
and that there are some oscillations, reflecting the shell structure as well, 
and inflexions points. Deviations from linearity, indicated by curvature, are 
more obvious when one investigates the relationship {D — u) versus u, which 
is illustrated in Fig. [H The amplitude of the oscillations in the {D, u) 
relationship is very small, which is not the case of (P, p) relationship. 

6 Conclusion 

Shock waves generating finite-temperature dense matter make possible the 
generation in laboratory experiments of extremely high energy densities typ- 
ical of matter in the few first microseconds after the creation of the universe 
and for such astrophysical objects as stars and giant planets. The physical 
information obtained from these experiments extends the basic knowledge 
of physical properties of these systems to pressures nine orders of magnitude 
higher than found, e.g., in our atmosphere. We proposed a qualitative and 
quantitative study of quantum orbital effects on the principal shock adia- 
bat for different elements, for two quantum self-consistent-field models: the 
atom in a spherical cell (ASC) and the atom in a jellium of charges (AJC). 
Quantum orbital effects lead to oscillations corresponding to the ionization 
of successive orbitals that are also visible in the electronic specific heat, and 
in the energies of the orbitals. The AJC model provides a better treatment 
of pressure ionization, since it relies on a full quantum treatment of the 
electrons. An estimate of the maximum compression, giving realistic values, 
has been proposed. 

The next step will be to test whether the oscillations still exist "beyond" 
the Average Atom model, i.e., when a variety of electronic configurations 
is taken into account [7j. It remains a difficulty in the existing models to 
represent, in a simple and suitable way, the infiuence of the plasma environ- 
ment on a specific ion. Indeed, as this environment fluctuates, the number, 
the localization in space and the structure of neighboring ions will change 
drastically. Therefore we will in the future represent the ionic environment 
including radial distribution functions, thereby going beyond the adiabatic 
approximation. 
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Figure 2: Electronic specific fieat for Al. 
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Figure 3: Maximum compression rate obtained from [22], from our model 
and compared to the maximum compression observed from a TF, ASC and 
AJC calculations for Be, B, C, Na, Mg, Al, Fe and Cu. 
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Figure 4: (D — u) versus particle velocity u for Be, Al and Fe. 
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